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Abstract. 

In the particle-flow approach information from all available sub-detector systems 
is combined to reconstruct all stable particles. The global event reconstruction 
has been shown to improve, in particular, the resolution of jet energy and missing 

(Jjj I transverse energy in pp collisions compared to purely calorimetric measurements. This 

improvement is achieved primarily by combining the precise momentum determination 

y , of charged hadrons in the silicon tracker with the associated energy depositions in the 

(-H ' calorimeters. By resolving individual particles inside jets, particle flow reduces the 

sensitivity of the jet energy scale to the jet fragmentation pattern, which is known to be 
one of the largest sources of systematic uncertainty in jet reconstruction. Particle flow 

^ ' reconstruction is thus potentially well-suited for the study of potential modifications to 

0^ , jet fragmentation in heavy-ion collisions. The particle flow algorithm has been adapted 

to the heavy-ion environment. The performance of jet reconstruction from particle flow 
objects in PbPb collisions using the anti-Zcx jet reconstruction algorithm is presented. 



The CMS detector [T] contains four primary sub-detector systems, which are, in 
order of increasing distance from the beam, a high precision silicon tracking system, 
a high resolution, high granularity electromagnetic calorimeter (ECAL), a somewhat 
/\ • coarser hadron calorimeter (HCAL), and finally muon chambers. Using a nearly 4 T 

c^ ■ magnetic field, this combination of sub-detectors, allows one, in principle, to distinguish 

five categories of stable or pseudo-stable particles: charged hadrons, photons, neutral 
hadrons, muons and electrons. 

In the approach traditionally used at hadron colliders, jets are reconstructed 
exclusively from calorimeter information (Calo jets). In CMS jets are clustered from 
calorimeter towers which consist of a single HCAL cell and the subtended 5x5 set of 
ECAL crystals. As the tower dimensions are determined by the HCAL geometry, the 
ECAL granularity is not exploited in this approach. The jet resolution is driven by 
the energy resolution of the HCAL {aE/E ^ 100%/ \/E). Moreover the particle px 
dependence of the calorimeter response is not corrected for, leading to a non-linear 
response to jets, causing the jet response to depend on the fragmentation pattern of 
jets, which translates into an increased uncertainty on the jet energy scale. To further 
complicate matters, at low p^, due to the strong magnetic field, charged hadrons are 
bent outside of a typical jet radius. 
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An alternate approach, known as particle flow (PF), individually reconstructs 
all stable particles by combining information from all available sub-detectors. These 
reconstructed particle may then be used as input to a jet reconstruction algorithm in 
place of calorimeter towers. By making use of the high precision silicon tracker, one may 
improve the p^ resolution of charged hadrons, and largely circumvent the non-linearity 
of the HCAL response. The CMS implementation of PF has been shown to largely 
improve the jet px and angular resolution in pp collisions [2]. In heavy-ion collisions 
the same PF algorithm is employed. However, charged particle tracks are reconstructed 
with a different algorithm to cope with the high multiplicity of PbPb collisions. A hit is 
required in each of the three layers of the pixel detector which keeps the fake rate very 
low, even for the most central collisions, but reduces the efficiency to a level of 60-70%, 
compared to about 90% for the pp algorithm. Charged hadrons for which no track is 
reconstructed default to a purely calorimetric measurement (corrected for the expected 
hadron response), degrading the jet energy resolution compared to the pp algorithm [2]. 

Jet reconstruction in PbPb collisions proceeds as follows. First, the underlying 
event is subtracted according to the procedure outlined in [3] and first employed in [3]. 
As is the case for pp collisions, the jet energy corrections ( JEC) are derived from pythia. 
Data-driven techniques in pp, such as 7+jet balancing, show that these corrections are 
valid to within a few %. The validity of this factorized background subtraction and JEC 
approach can then be tested in pythia events embedded in PbPb events simulated with 
HYDJET (pythia-|-hydjet) |5]. Jets are reconstructed with the anti-fcx algorithm with 
R = 0.3. Figure fTal shows the jet reconstruction efficiency in 10% most central hydjet 
events for both PF jets and Calo jets reconstructed with a cone algorithm with R = 0.5 
as described in [4]. The PF jet reconstruction is fully efficient by about 50 GeV/c, a 
somewhat lower value than for Calo jets, mostly owing to the smaller value of R. 
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Figure 1: (a) The jet reconstruction efficiency in 10% most central HYDJET events. 



(b) Unmatched jets per leading jet in 10% most central pythia -(-hydjet events. 
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For an analysis of dijets, it is important not only to find jets, but to be able to select 
jet pairs from the same hard scattering. This issue is investigated in pythia+hydjet. 
First, a leading jet is selected with px > 100 GeV/c. We then count the rate at 
which reconstructed jets on the away-side cannot be matched in AR to an embedded 
PYTHIA jet. Figure flbl shows the rate of unmatched jets per leading jet as a function of 
reconstructed jet pt- The PF jet selections shows a very small rate of unmatched jets at 
40 GeV/c, significantly lower than Calo jets at the same px, which is at least partially 
driven by the smaller value of R. 
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Figure 2: The inclusive jet response and resolution for PF and Calo jets for raw 
(uncorrected - open symbols) and corrected (closed symbols) jets. 



Figure 2 further demonstrates the performance of PF jets in pythia+hydjet, 
comparing the reconstructed jet px to that obtained by clustering generator level 
particles. The top row shows the jet px resolution, extracted using a Gaussian fit, for 
several different centrality selections. PF reconstruction improves the jet resolution over 
most of the px range. This improvement is, however, partially offset by the decreased 
value of R, which increases the degree to which energy can migrate into and out of 
the jet, when comparing the reconstructed and generator-level jets. The bottom row 
of Fig. [2] shows the jet response, defined as the mean of the ratio of the reconstructed 
to generated jet px, both before and after the application of JECs. As expected, the 
raw response to PF jets is closer unity, showing that the size of the required JECs is 
reduced by about half. The corrected response nicely closes at unity independent of 
centrality. This demonstrates the validity of the combined background subtraction + 
JEC procedure and shows that the particle-based jet reconstruction is insensitive to 
occupancy effects. 
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As the JECs corrections are based on an inclusive jet selection in pythia, the 
true response may differ from that obtained from simulation. This may arise if the 
fragmentation of jets is not accurately modeled, for example if the fraction of quark 
and gluon jets is different in simulation and data. This may be of particular concern in 
heavy-ion collisions where quenching is expected to be larger for gluons than for quarks, 
and the parton distribution functions are less well known than for protons. In Fig 3 
the corrected response is shown differentially for quark and gluon jets. One observes a 
splitting between the quark jets, which tend to fragment harder than average, and gluon 
jets, which fragment softer. However, this splitting smaller for PF jets than Calo jets, 
indicating that PF jet reconstruction is less sensitive to the details of the fragmentation. 
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Figure 3: The corrected jet response for quark and gluon jets with PF and Calo jet 
reconstruction in PYTHIA. 



To conclude, PF reconstruction has been performed in heavy-ion collisions for the 
first time. The jet finding efficiency and the dijet mismatch rate have been evaluated in 
simulated data. The performance of PF-based jet reconstruction has been demonstrated 
both in terms of both response and resolution and compared to Calo-based jet finding. 
Finally, by studying the response of quark and gluon jets, PF jets were shown to be less 
sensitive to the details of jet fragmentation than Calo jets. 
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